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ABSTRACT 

The best measured rotation curve for any galaxy is that of the dwarf 
spiral DDO 154, which extends out to about 20 disk scale lengths. It provides 
an ideal laboratory for testing the universal density profile prediction from 
high resolution numerical simulations of hierarchical clustering in cold dark 
matter- dominated cosmological models. We find that the observed rotation 
curve cannot be fit either at small radii, as previously noted, or at large radii. 
We advocate a resolution of this dilemma by postulating the existence of a dark 
spheroid of baryons amounting to several times the mass of the observed disk 
component and comparable to that of the cold dark matter halo component. 
Such an additional mass component provides an excellent fit to the rotation 
curve provided that the outer halo is still cold dark matter-dominated with 
a density profile and mass-radius scaling relation as predicted by standard 
CDM-dominated models. The universal existence of such dark baryonic 
spheroidal components provides a natural explanation of the universal rotation 
curves observed in spiral galaxies, may have a similar origin and composition to 
the local counterpart that has been detected as MACHOs in our own galactic 
halo via gravitational microlensing, and is consistent with, and even motivated 
by, primordial nucleosynthesis estimates of the baryon fraction. 

Subject headings: dark matter — galaxies: formation and halos 
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1. Introduction 

Cosmological cold dark matter theories of structure formation in the universe via 
hierarchical merging are in some difficulty. Recent high-resolution cosmological N-body 
simulations (Navarro et al. 1996a, 1997) have substantially improved our understanding 
of the equilibrium density profiles which dark matter halos achieve when formed through 
hierarchical clustering. It has been shown that the violent, coUisionless dynamical relaxation 
processes during the formation phase of dark matter halos lead to equilibrium profiles that 
have similar shapes, independent of halo mass, initial density fluctuation spectrum, and 
adopted cosmological model. All dark matter profiles can be well fit by the simple formulae 
for density and mass: 

MoMir) = Mox[ln{l + r/R,)-{r/Rs)/{l + r/Rs)], (1) 

where po is the density of the dark matter halo evaluated at the scale radius Rs and Mq 
is the characteristic mass; Rg and Mq (or po) ^-re free parameters. Mq is a function of the 
total virial mass M200 inside the virial radius i?2oo 

M = (2) 

° ln(l + R200/Rs) - {R200/Rs)/{1 + R200/Rs) ^ ' 

where i?2oo denotes the radius inside which the averaged overdensity of dark matter is 200 
times the critical density of the universe. The simulations also show that, for any particular 
cosmology, Rg and Mq are strongly correlated: Rg = 1.63 x lO~^~^(M2oo/-^0)^^^^~^^^ kpc 
where c = log(it!2oo/-Rs) is approximately 1.4 for low-mass dark matter halos as considered 
in this paper. Low-mass halos are denser than more massive systems. This results from the 
fact that lower mass halos form earlier, at times when the universe is significantly denser. 
Dark matter halos therefore represent a 1-parameter family, being completely described 
by equation 1 and their virial mass M200 or virial radius i?2oo- The universal profile has 
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been verified by simulations to halo masses as small as M200 ~ IO^^Mq but there is no 
reason to believe that these results would not be valid for halos which are even one order of 
magnitude lower in mass. 

Substantial progress has also been made within the past decade on the observational 
front. The observations show that the rotation curves of low-luminosity disk galaxies and 
low-surface brightness galaxies are strongly dark matter-dominated (Persic & Salucci 1988, 
1990, Puche & Carignan 1991, Broeils 1992, Persic et al. 1996, de Blok et al. 1996). 
The spherically-averaged mass distributions M(r) = V^r/G can be derived from the 
measured circular velocity distributions V^(r) of gaseous galactic disks. Subtraction of the 
contributions by the visible components gives the dark matter mass profiles Mdm{t) or 
the corresponding dark matter rotation curves Vc,DM{r) = {GMdm/'^Y^'^- Given these 
theoretical and observational developments, an important question arises as to whether the 
theoretically predicted rotation curves, determined by equation 1, are in agreement with 
the observations. 



2. The dark matter halo of DDO 154 

The low-surface brightness dwarf galaxy DDO 154 has one of the most extended and 
best-studied dark matter halo rotation curves (Carignan & Freeman 1988, Carignan & 
Beauheu 1989, Carignan & Purton 1997), with a precise decomposition into contributions 
from stars, gas and dark matter. It is also one of the most gas-rich systems known with an 
inner stellar disk component of mass 5(±2.5) x l(f Mq and an extended HI disk with scale 
radius of 0.4(±0.05) kpc and mass of approximately 3(±1) x lO^M© (Carignan & Beaulieu 
1989). The shape of the rotation curve, even in the innermost regions, is completely 
dominated by the dark matter halo with a total mass within i?2oo of about 5 x lO^M© (see 
section 3). 
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DDO 154 is a good candidate with which to test cosmological models. Fig. 1 (upper 
panel) shows the dark matter rotation curve of DDO 154 and compares it with the profiles 
predicted from cosmological models. The stellar disk and the HI contribution (Carignan 
& Freeman 1988) have been subtracted, adopting a stellar mass-to-light ratio of (M/Lb)* 
= 1. The error bars indicate the uncertainties in the observations. Note that the data 
extends out to 21 disk scale lengths and that the rotation curve clearly decreases beyond 5 
kpc. The thick dashed and dotted curves fit the theoretically predicted rotation curves as 
determined by equation 1 to the innermost and outermost regions, respectively. Fitting the 
inner regions (dotted curve, labeled A) is well known to pose a problem (Flores & Primack 
1994, Moore 1994, Burkert 1995) as the theoretical models predict a central density 
cusp, whereas the observed velocity profile indicates a large isothermal core with a constant 
density. As a result, the theoretical models lead to far more mass in the innermost region 
than is seen. It has been suggested that this discrepancy could be solved by assuming 
secular processes (e.g. violent galactic winds) in the baryonic component (Navarro et al. 
1996b) which could also affect the innermost parts of dark matter halos. 

However a similar problem exists in the outermost regions. Whereas the observed dark 
matter rotation curve clearly decreases beyond 5 kpc, the theoretically predicted universal 
profile fit leads to a very massive and extended halo (M200 = 1-8 x IO^^H'^Mq; -R200 = S63h~^ 
kpc, where h is the Hubble constant in units of 100 km s~^ Mpc~^) with a rotation curve 
that increases beyond 8 kpc. The outer regions certainly cannot be affected by secular 
mass loss involving a baryonic component, and a different explanation must be sought. 
The dashed line in the upper panel of Fig. 1 (labeled B) shows a fit to the outer rotation 
curve of DDO 154. In this case the dark matter excess in the inner regions is unacceptably 
large.Approximately lO^M© in dark matter would have to be moved from the inner 2 kpc 
into the region between 2 kpc and 4 kpc to explain the discrepancy between the observed 
inner rotation curve and the predicted one. This is far more than would be expected as a 
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result of secular processes. 

Figure 2 compares the values of Mq and Rs (see equation 1), derived from fitting the 
inner or outer rotation curve of DDO 154 with the standard, cluster-normalized CDM 
predictions (Navarro et al 1996a, 1997). Note that both fits require pairs of parameter 
values which are not in agreement with theory. In order to lie within the expected parameter 
range one would have to choose values for Mq and Rg which do not fit the rotation curve, 
either in the inner or in the outer regime. 

3. The dark bciryonic component of DDO 154 

What is wrong with either CDM or with the failure to fit the rotation curve of DDO 
154? There is accumulating evidence that the CDM models can account for many aspects 
of large-scale structure and therefore should not be dismissed. Errors in determining the 
rotation curve of DDO 154 are equally unlikely due to its very regular velocity field from 
which Vc{r) can be derived and unambiguously corrected for warping. As described above, 
secular processes cannot explain the discrepancy. There also cannot exist much more mass 
in the HI disk, e.g. in H2, as such a massive disk would be gravitationally unstable and 
would efficiently form molecular clouds and stars, in contradiction to the observations. 

We propose, that in addition to the CDM dark matter halo which can be represented 
by the one-parameter family of universal profiles, DDO 154 contains a second dark but 
compact baryonic component whose presence we infer in order to explain the observed 
rotation curve. Because of the arguments mentioned above, this new component cannot be 
gaseous and located in the disk. We identify it with a spheroidal distribution of massive 
compact baryonic objects, MACHOS, located in the inner regions of the dark matter halo. 

There is in fact ample room for such a subdominant baryonic component, relative to 
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the massive dark halos of the standard CDM theory in the form of MACHOS. Primordial 
nucleosynthesis requires a baryonic component of f« 0.015 ± 0.008 (Kurki-Suonio, 

Jedamzik & Mathews 1997, Copi, Schramm & Turner 1995) whereas the observed value 
for stellar and gaseous components in disks hes in the range of Qd ~ 0.004 ± 0.002 (Persic 
& Salucci 1992, note that Qd ~ 0.006 for DDO 154, assuming Qhaio = 0.1). Where are all 
these additional baryons? Most of them should be found in galaxies unless winds have 
strongly depleted all the initial baryons during the protogalaxy phase. Even if protogalactic 
winds have occurred, as is suggested by theory and observation, it is unlikely that all of the 
baryons that are not in the disk would have been lost. The recent MACHO experiments 
(Alcock et al. 1996) demonstrate that there indeed exists a substantial baryonic component 
in a more extended, spheroidal distribution in the Milky Way, at least between the Sun and 
the LMC, in addition to the observed stellar and gaseous baryonic components. It is hkely 
that a similar, yet hitherto unobserved, component also exists in DDO 154. 

The mass distribution Msphir) of the proposed spheroidal dark baryonic component is 

determined by 

MsAr)-V,'^^^r/G-M*^^{r) (3) 

which is the difference in total mass inside a radius r as expected from the dark matter 
rotation curve Vc^dm and the dark matter mass M^j^ inside r. We consider two extreme 
possibilities for determining M^j^{r). One possibility is just to adopt the universal profile 
predicted by hierarchical infall models (equation 1): M^^(r) = M^Mi^)- In this case 
the baryonic dark component formed at the same time or even earlier than the extended 
dark halo. In the opposite limit we imagine that the MACHO spheroid formed after the 
extended dark halo profile was established. In this case, if the dark baryonic component is 
of comparable mass inside a certain radius, the dark halo within this radius will undergo 
adiabatic contraction. For simplicity we will assume spherical symmetry and a formation 
timescale of the MACHO halo which is long compared to the dynamical timescale of the 
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dark halo. In this case the corrected dark matter mass profile is determined by (Binney & 
Tremaine 1987) M^^(r) = Mouir*) where MnMir) is determined by equation 1 and 

MDM{r*) ■ ^ ^ 

Inserting Msph (equation 3) into equation 4 and choosing values for Mq and Rg, one can 
determine iteratively the adiabatically contracted dark matter mass profile M^j^^ and 

Our approach seems at first sight a very contrived solution to the problem of the 
rotation curve, as we attribute the discrepancy between observations and theory to a third, 
invisible component which introduces an extra degree of freedom. How can this idea be 
tested? There exist in fact very strong constraints on the shape of Msph{r), which must 
be fulfilled in order for our model to be physically correct. First of all, the total mass 
distribution of the dark spheroid has to be positive everywhere and increase monotonically 
with increasing radius, up to a maximum radius Rgph beyond which it has to remain 
constant. Second, the density of the spheroid Psph{r) = {dMsph/ dr) / {ATir"^) must decrease 
monotonically with increasing radius. 

These constraints require that the dark matter rotation curve lies everywhere below 
the observed rotation curve. Dark matter fits like the curves labeled A and B in the upper 
panel of Fig. 1 are ruled out. Rgpu could in principle lie outside the observed radius regime. 
However for DDO 154 the measured rotation curve extends out to 21 disk scale lengths. 
We identify the dark spheroid with an intermediate baryonic component that presumably 
formed during the dissipative protogalactic collapse phase, and which should therefore be 
more centrally concentrated than the dark matter halo. The rotation curve of DDO 154 
clearly decreases outside 5 kpc, and this indicates that a large fraction of its mass lies at 
smaller radii. As the dark matter halo dominates the total mass in the outer regions, we 
expect that the mass distribution of the baryonic spheroid becomes constant and equal to 
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its total mass M*"^ at a radius Rgph < 5 kpc. 

Note that the radial dependence of Vc and Mdm is given by the observations and 
equation 1. It is therefore not trivial that the difference of the terms on the right hand side 
of equation 3 gives a profile which meets all of the constraints mentioned above, even if the 
dark matter scale parameters Mq and Rg are treated as free parameters. We have indeed 
found such a solution, however only for very special values of Mq and Rs, values which for 
the models with and without adiabatic contraction lie within the small shaded and dark 
regions of figure 2, respectively. These sets of solutions uniquely determine the total mass of 
the spheroid M^^f^, and hence, via equation 2, the total mass M200 of the dark matter halo. 

The solid lines in the upper and lower panels of figure 1 provide an excellent fit to 
the rotation curve, adopting our Macho model with the non-baryonic dark matter halo 
parameters fixed at Mq — Ax 10^ Mq and Rs — 4.5 kpc for the case without adiabatic 
contraction and Mq = 3 x lO^M© and Rg = 5.5 kpc if adiabatic contraction is included. The 
inserts within Figure 2 show the mass and density profiles of the dark baryonic spheroid 
and the dark matter halo without and with adiabatic contraction. Note that in the former 
case without adiabatic contraction, the MACHO component would dominate the inner mass 
and density distribution, whereas in the latter case the baryonic and non-baryonic halos 
have comparable masses in the inner region. For all acceptable models we find Rgph ~ 5 
kpc, Ml^l f« 1.5 X 10^ Mq and M200 ~ 5 x 10^ Mq. The density profiles of the MACHO 
component can be well approximated by an isothermal sphere with a constant velocity 
dispersion of cr = 40 km/s. 

There exists an additional independent constraint in order for our model to be 
acceptable within the framework of standard cosmology: Mq and Rg must follow the tight 
relationship, predicted by cosmological models. Indeed, figure 2 shows that the values 
which have to be adopted in order to give a physically correct mass distribution overlap 
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with the parameter space of values expected from standard cosmology, in contrast to 
the one-component dark matter fits (starred points). Moreover, as we identify the dark 
spheroid with the missing baryonic component, early universe nucleosynthesis requires that 
n,/na = (M*^^ + Mlt,,i:)/Mlt^,i^ ^ 5(±3) x h-\ where M*^*,,,, = 3 x IO^Mq is the total 
visible mass of DDO 154. For the areas of possible solutions we find a dark-to-luminous 
mass fraction of approximately 6, in agreement with the expectations. 

Note that an arbitrary choice of Mq, Rg and Mgp^ would almost certainly fail to 
meet all of these constraints. The excellent agreement of our model with the predictions 
from cosmological models of structure formation and primordial nucleosynthesis provides 
additional evidence for the presence of a dark baryonic component in DDO 154. 

4. Discussion and Conclusions 

Our three component model of luminous baryons in a disk configuration, and 
MACHOS and cold dark matter in a spheroidal distribution, can reconcile the most detailed 
observations of a rotation curve to date with the hierarchical clustering theory of galaxy 
formation. This might be the first (indirect) detection of a MACHO component in another 
galaxy. It also has allowed us to study in detail for the first time the internal density 
distribution of a dark baryonic spheroid due to the excellent high-resolution data of DDO 
154's rotation curve. 

The structure of the MACHO spheroid in DDO 154 is surprisingly similar to the 
MACHO halo of the Milky Way, the existence and mass of which has been inferred from 
a completely different method: gravitational microlensing events of stars in the Large 
Magellanic Cloud. Alcock et al. (1996, 1997) find for the Galaxy that within 50 kpc (14 
disk scale lengths) the total masses of MACHOS and dark matter are comparable and of 
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order 2.5 x lO^^Mg. This is 4 to 5 times the mass of the galactic disk (M^ f« 6 x 10^°Mq). 
The dark baryonic spheroid of DDO 154 also extends out to 14 disk scale lengths at which 
radius the mass of the MACHO halo is again similar to the mass of the dark matter halo, 
with, in this case, a mass of order 1.5 x lO^M©. The inferred MACHO mass is also of order 
5 times the mass of the HI disk. This agreement in the relative mixture of dark matter, 
MACHOS and disk material indicates that these components formed in both galaxies 
from a similar continuous dynamical process, with the MACHO spheroid representing a 
presumably dissipative component intermediate between the coUisionless non-baryonic dark 
halo and the strongly coUisional, dissipation-dominated, rotationally-supported disk. This 
might provide an explanation for the puzzling observational result that disk galaxies have 
universal rotation curves (Casertano & van Gorkom 1991, Rubin et al. 1985, Persic et 
al. 1996), requiring a connection between their galactic disks and their dark spheroidal 
components. Universal rotation curves would be expected for galaxies of any mass where 
the relative mass and radius ratios between the dark matter halo, the MACHO halo and 
the disk are universal numbers. 

Our results indicate that the baryonic component of DD0154 and probably also of 
other disk galaxies consists of two components, a spheroidal MACHO component which 
represents about 22% of the total mass and a visible disk component with only 4% of the 
total mass. The total baryon fraction in galaxies is then of the order of a quarter of the 
total mass, a value which is higher than expected from the primordial nucleosynthesis 
predictions if the relative mixture of baryonic and non-baryonic matter is universal. This 
baryon segregation could result from dissipative processes during the formation of halos 
which concentrate the baryons relative to the nondissipative dark matter. 

The origin of two separate baryonic components, namely a dominant dark spheroidal 
component and a disk component, is an interesting and yet unsolved theoretical puzzle 
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which could provide important information on the dissipative formation history and 
evolution of galaxies. 
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FIGURE CAPTIONS 

FIG. 1. Upper panel: The dark matter rotation curve of DDO 154 is shown with error 
bars. The dotted hne (A) shows a fit to the inner parts of the rotation curve, adopting the 
dark matter halo structure as predicted by cosmological models. The dashed line (B) shows 
a dark matter halo fit to the outer part of the rotation curve. The solid hne shows the fit 
achieved with the 2-component MACHO model without adiabatic contraction, assuming 
dark matter halo parameters Mq = 4 x lO^M© and Rg = 4.5 kpc. The lower dashed and 
dot-dashed curves show the contribution to the rotation curve of the MACHO spheroid and 
the dark matter halo (C), respectively. 

Lower panel: The 2-component MACHO model with adiabatic contraction is shown, 
adopting halo parameters Mq = 3 x IO^Mq and Rg — 5.5 kpc. The dashed curve shows 
the contribution by the MACHO spheroid. The dot-dashed curve (D) and the dotted curve 
(C*) show the contribution of the dark matter halo after and before adiabatic contraction, 
respectively. 

FIG. 2. — Standard cluster normahzed cold dark matter models predict that the dark 
matter halo scale radii Rg and scale masses Mq should lie within the narrow band enclosed by 
the two parallel sohd lines. The parallel dashed hues enclose the region of scale parameters, 
expected for the less favoured COBE-normalized cold dark matter model. One-component 
dark matter fits to the rotation curve of DDO 154 would result in scale parameters as shown 
by the two stars for the inner (A) and outer (B) fits. The two-component MACHO model 
without adiabatic contraction and with adiabatic contraction requires the scale parameters 
to lie within the dark area (labeled C) and the shaded area (labeled C*), respectively. 
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The upper inserts show the mass and density distribution (p in units of Mq pc~^) for the 
standard model (Mq = 4 x 10^Mq,Rs = 4.5 kpc, star inside dark area) without adiabatic 
contraction, with the sohd and dashed hnes representing the dark baryonic spheroid and the 
dark matter halo, respectively. The lower insert shows the mass distribution of the standard 
model (Mq = 3 x IO^Mq, Rg — 5.5 kpc, star inside shaded area) with adiabatic contraction. 
The dots represent the total dark matter mass profile as predicted from the rotation curve. 
The lower sohd line shows the mass profile of the MACHO halo. The dot-dashed and the 
dotted curves show the mass distribution of the non-baryonic dark matter halo after and 
before adiabatic contraction, respectively. 



